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Abstract 

O |. We study the limiting thermodynamic behavior of the normalized sums of spins in multi-species 

^ r*| Curie-Weiss models. We find sufficient conditions for the limiting random vector to be Gaussian (or 

-I— > . 

to have an exponential distribution of higher order) and compute the covariance matrix in terms of 
model parameters. 

(N 
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Introduction 

The study of the normalized sum of random variables and its asymptotic behavior has been and contin- 
ues to be a central chapter in probability and statistical mechanics. When those variables are independent 
and have finite variance the central limit theorem ensures that the sum with square-root normalization 
converges toward a Gaussian distribution. The generalization of that result to non-independent vari- 
ables is particularly interesting in statistical mechanics where the random variables have an Hamiltonian 
interaction. 

Ellis and Newman [IH3] have studied the distribution of the normalized sums of spins whose interaction 
is described by a wide class of mean field Hamiltonian a la Curie- Weiss. They have found the conditions, 
in terms of the interaction, that lead in the thermodynamic limit to a Gaussian behavior and those which 
lead to a higher order exponential probability distribution. 

In recent times a multi-species extension of the Curie- Weiss model has been proposed in the attempt 
to describe the large scale behavior of some socio-economic systems [4] . Multi-populated non- interacting 
spin models are at the basis of the so called Mc Fadden discrete choice [5] theory. The extension of the 
discrete choice theory to the interacting, and more realistic, case is an important problem toward the 
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understanding of the collective behavior of social and economical systems. The investigation of multi- 
species models has been pursued at a mathematical level in [Bj where there have been proved properties 
like the existence of the thermodynamic limit by monotonicity, the computation of the free energy and 
of the intensive quantities like local magnetizations. The phenomenological test of the model has been 
started in [7] and it is a topic of current investigations. 

In this paper we deal with the study of the normalized sum behavior for a multi-populated model 
with mean field interaction. We prove that, under the assumption that the mean field Hamiltonian 
interaction has a convexity property, when the system reaches its thermodynamic limit the random 
vector whose components are the sums of spins on each population, converges to a nontrivial random 
variable S. The behavior of S depends crucially upon the nature of the minima points of a function 
G (the pressure functional) which we associate to the model interaction type. In particular it is the 
value of the determinant of the Hessian matrix of G computed on the minima points that establishes 
the Gaussian or non Gaussian behavior of the random vector. In the case of a unique minimum point 
if the determinant is different from zero then S is a multivariate Gaussian whose covariance that can 
be computed from the mean field equations. Otherwise S has a distribution whose density behaves as a 
higher order exponential. When the function G has more than one minimum point we obtain a similar 
result whenever the random vector of the magnetizations is close enough to one of the minimum points. 

This work is organized as follows. Chapter one introduces the language and the notations and states 
the main results in theorems Q] and [2] Chapter 2 contains the proofs. Chapter 3 describes specific cases 
in which the distribution is Gaussian and others in which is not. The appendix contains the proof of the 
lemmas that make the paper self contained. 

1 Definitions and Statements 

We consider a system of N particles that can be divided into n subsets P±, . . . , P n with Pi n P s = 0, 
for I 7^ s and sizes \Pi\ = JVj, where YaLi W = N. Particles interact with each other and with an external 
field according to the mean field Hamiltonian: 



The Oi represents the spin of the particle i, <Ji = ±1 while Jij is the parameter that governs the mutual 
interaction between the particle i and the particle j and takes values according to the following symmetric 
matrix: 
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where each block 3i s has constant elements J; s . For Z = s, 3u is a square matrix, whereas the matrix 
3i s is rectangular. We assume Jn, J22, • ■ • , Jnn be positive, whereas J; s with Z ^ s can be either positive 
or negative allowing both ferromagnetic and antiferromagnetic interactions. The vector field takes also 
different values depending on the subset the particles belong to: 



ATi{ 
N 2 { 



h 2 



v / 



where each h; is a vector of constant elements hi . 

The joint distribution of a spin configuration er = (01, . . . , ctjv) is given by the Boltzmann-Gibbs 
measure: 

1 N 

Pn,jA°} = ^rnr^ HN(rT) II d M ( 2 ) 

where Z N (3, h) is the partition function 

N 

Z N (3,h) = J e- H "^]ldp(<Ti) 

and p is the measure: 



»=i 



(3) 



(4) 



where S(x — x ) x € R denotes the unit point mass with support at x . The inverse temperature j3 isn't 
explicitly written because it is included in the parameters of the Hamiltonian. 
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By introducing the magnetization of a set of spins A as: 



m A 



1 1 lEA 



(5) 



and indicating by to/(ct) the magnetization of the set Pi, and by ai — Ni/N the relative size of the set 
Pi, we may easily express the Hamiltonian ([!]) as: 



where the function g is: 

g(mi(er), m n (<x)) = - ^ aia s Ji s mi{o-)m s (cr) + ^ aihimi(tr) . 

1,8=1 1 = 1 

In [6] it is shown that the thermodynamic limit of the pressure function 

r(J,h) = iln($>-*M^ 



PN\ 



(6) 



(7) 



(8) 



exists and is reached monotonically if the the function g given by (J7]) is convex (see also (§]). In this case: 



lim pjy(J,h)= sup p(x) 



JV— >oo 



(9) 



where the functional p is: 

^ n n / / n 

) = In 2 - - ^ aia s Ji s xix s + a x In I cosh I ^ a s J ls x s + hi ) ) . (10) 



1,8=1 1 = 1 ' X S=1 

It is easy to verify that the function g is convex if the following matrix, called reduced interaction 
matrix, 



Jn J\2 

J12 J22 



Jln^ 
J211 



(11) 



\Jln J271 • • ■ Jnn J 

is positive dchnitc. 

The extremality conditions of p give the Mean Field Equations of the model 



pi = tanh I a iJu M + h 
v 1=1 

fi 2 = tanh ( aiJ 2 i pi + h 2 ) 
v 1=1 J 



(12) 



fj, n = tanh 



V 1=1 



In the thermodynamic limit the random vector (mi(<r), . . . , m ra (<r)) weakly converges, with respect 
to the Boltzmann-Gibbs measure, to the deterministic vector (pi, . . . , fj, n ) solution of the Mean Field 
Equations. This means that the variances of the magnetizations vanish for large N a part on isolated 
singularities (see [5] for the precise statement). In this paper, define the sum of the spins of a set A as: 

S A (<r) = J2 a i ( 13 ) 

and indicating by Sz(er) the sum of the spins of the set P/ we want to determine a suitable normalization 
for the vector elements Si(<x), . . . , S n (cr) so that in the thermodynamic limit they converge to well defined 
random variables with finite (non zero) covariance matrix. The problem in n = 1 has been solved in pQ 
and [3]. 

We shall see that the behavior of the limiting distribution of the sums of spins depends crucially on the 
number and the nature of the minima points of the function G = —p + In 2 



G(xi, ...,£„) = - ^ aia s Ji s xix s - ^ «; In ( cosh ( a s Ji s x s + hAj . (14) 

l,s=l 1=1 ^ ^s=l ' ' 

Let fj, 1 , . . . , fj, p be global minima points of the function G. For each p there exist the functions 



(x) > 0, such that around fi p we can write G as: 



d / j\ 

G(x) = Gin?) + ]T G% (x - v?) + o(J|x' - n*>'\ 2 + |x" - /xp"| 2 ^) J 



where (x',x") is a partition of the coordinate x and q is a positive rational number such that 1/q € N 
and 

Gf (ix', i 9 x") = t 2j G^ (x', x") all t > 0. (16) 

We define the type k p of the minimum point fi p as the smallest d such that G% d (x — /x p ) ^ as x ^ 
and Gj) (x — /x p ) = for j = 1, . . . , d — 1. We observe that when q = 1 the expression (JTSj is the 
Taylor expansion of the function G. In this case, the only one we treat in this paper, k p is called the 
homogeneous type of the minimum point fi p . In particular if a minimum points fi p has homogeneous 
type equal to 1, around y, we have: 

G(x) = G(n p ) + ^JC G (M p )(x - » p ), (x - H p )) + o(||(x - ^) 2 ||) (17) 

where IKg(/.* p ) is the Hessian matrix of G computed in the minimum point fi p . 
We introduce some useful notations. Considering x, y £ K™ and 7 G R we define 

• x 7 = (x 7 , . . . , x 7 ); 

• xy = (xij/i, . . .,x n y n ); 
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(-L 



where yi 7^ for 1 = 1, 



x n 

y ^yi y n ■ 

Now we can state our main results. The following theorem solves the problem of the correct normal- 
ization of the random vector of the sums of spins whenever the function G admits a unique minimum 
point. 

Theorem 1. Consider the mean-field Hamiltonian Hn = —N g{rn\(cr), . . . , m„(er)) where g is the convex 
function defined in Let fi = (fj,i, . . . , /i n ) be the unique global minimum point of the function G given 
by J_?^[ ). Let k be the homogeneous type of the minimum point. 



1. If k = 1 the random vector 



Si(<r)-Nifii S n (a)-N n n r , 



(18) 



as Ni — > 00, . . . , N n — > 00, for fixed values of a±, . . . ,a n , weakly converges to a normal multivariate 
distribution whose covariance matrix is: 



( 



X = 



dm 

dhi 



dh 2 dhi 



dhi dhi 



dh 2 



dh n dhi 
dh n dh 2 



dfii dfi n 
\\J dh n dhi 



d[i 2 d[i n 
dh n dh 2 



d[i n 
dh n 



(19) 



2. If k > 1 the random vector 

S k {<r 



Si{cr) - Ni^i S n (cr) - N n ^ r 



(20) 



as Ni — > 00, N n — > 00, for fixed values of ai, a n , weakly converges to a distribution 
proportional to: 

exp (-^-(^k)) ( 21 ) 

where a = (ai, . . . , a n ). 

The following second theorem handles the case in which the function G reaches the minimum in more 
than one point. 

Theorem 2. Consider the mean-field Hamiltonian Hn = —Ng(mi(cr), . . . , m„(er)) where g is the convex 
function defined in Let fi = (/ii, . . . ,/i n ) be a global minimum point of the function G given from 
Let k be the homogeneous type of the minimum point. Define S to be the minimum distance between 



G 



all distinct pair of global minimum points of the function G. Then for any d € (0, S) when the random 
vector of the magnetizations (mi(<x), . . . , m n {a)) is inside the ball B([jb, d) centered in /x of radius d 

1. if k = 1 the random vector ^(er) defined in U8\) as N\ — > oo, . . . , N n — > oo, for fixed values of ct\, 

a n , weakly converges to a normal multivariate distribution whose covariance matrix is given 

by W; 

2. if k > 1 the random vector Sk(cr) defined in A20\) as N\ — » oo, . . . , N n — > oo, for fixed values of ct\, 
. . . , a n , weakly converges to a distribution proportional to: 

-P(- G ^(^k)) P2) 

where a — (an, . . . ,a n ). 

2 Proofs 

2.1 Proof of Theorem [I] 

To prove the theorem we need the following lemmas. 

Lemma 1. Suppose that for each N , = (X^ , . . . , X^ ) and yjv = (Y^\ . . . , Y^ 1 ) are independent 
random vectors. Suppose that x^r weakly converges to a distribution v such that 

J e J<r ' x W(x) f for all reM". (23) 

R™ 

Then yjv weakly converges to /i if and only i/xjv + Yn weakly converges to the convolution is * fi of the 
distributions v and \x. 

This result follows because the weak convergence of measures is equivalent to pointwise convergence 
of characteristic functions (see |10j). 

Lemma 2. // the reduced interaction matrix J of a model defined by the Hamiltonian (0|) is positive 
definite, then the function G given in {1$ has a finite number ( different from zero ) of global minimum 
points and for any N 6 {1, 2, . . . } 

J exp ^— NG(x%, . . . ,x n )^jdxi . . .dx n < oo (24) 

R" 

See appendix for the proof. 

Lemma 3. Let A = D a JD a be a positive definite matrix where the matrix D a — diag{y/ai, . . . , y/a n } 
and the matrix J is defined in ill]) . Given the random vector (W\, . . . ,W n ) whose joint distribution is 
the normal multivariate 

Pi*) = ex P (- ^( Ax > x >) (25) 

7 



if (Wx, . . . , W n ) is independent of (Si(cr), . . . , S n (tr)) then for (fix, . . . , fx„) € 1" and 7 £ 1 the joint 
distribution of 



Wx W n \ (Sx{<t)-Nx»x S n (a)-N nPn 



is given by 



(26) 



cxp ( - NG^-p^ + nx, . . . , + Unjjdxx ■ ..dx n 
exp f- NG^-^y + fix,..., -^py + fir^jdxi . ..dx n 



(27) 



where G is the function defined in 
See appendix for the proof. 

We remark that as 7< 1/2, the random vector (Wx, ■ ■ ■ , W n ) does not contribute to the limit of (|2"7|) 
as Nx — > 00, . . . , N n — > 00. 

To prove the previous theorem we proceed as follows. For k > 1, by lemmas [T] and [3] with 7 = l/2k, 
we have to prove that, for any bounded continuous function ijj{x) : M n — > R 



L exp (- " g (n5s + ^))^ (x)dx 1 exp ( - g - (^w 



tp(x)dx 

(28) 



L 6XP (" iVG ( + #*) ) £ exp ( - (^) ) rfx 
where to easy the notation we set N = (Nx, . ■ ■ , N n ). Defined the function: 

B(x;/i) = G(x+/i)-G(/i) (29) 
there exists S > sufficiently small so that, as Nx — > 00, . . . , iV„ — > 00 for ||x/N 1 / 2 ' c || < (5 



x 



N ■ B ( n5» > " J = G 2, ( ) + W 



(30) 



+ o(l)P 2fe -i(x) 



where P2^(x) is a polynomial of 2fc order and P2^_i(x) is a polynomial of 2k — 1 order. 

Defined / = G(fi), for any closed subset V of M. n which does not contain fx there exist e > 0, so that as 

N -> 00 

e»f [ e- Ifa Mdx = o(e- If <). (31) 
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We pick S > as in ([30]) ■ By ([21]) there exists e > so that 

^ exp + ^^((II^) e ~"') ^ 



whereas by (|30p and dominate convergence, we have that 



/^^-P (- NG^ + = e-(/-(M)) /^^exp (- ^(n^)) ^(x)dx 

-/ R ex P (-G-(-^))^(x)dx. 



This proves the statement (2) of the theorem. 
We observe that for fc = 1 



^ fc (^)=^/(M)x,x) 



where (Kg = D a 1 JfcD a 1 is a positive definite matrix. In analogous way, we prove that for any bounded 
continuous function ^>(x) : M™ — > R: 



ex P ( - ^(^g(a*)x,x)) ip(x)dx. (33) 

The multivariate Gaussian obtained is the convolution of the distribution of the random vector (W\ , . . . , W n ) 
with the distribution of the random vector S 1 (cr) given by (|18|) . Indicated with <^>w(A), <^>gi(A) and 
w+ gi (A) respectively the characteristic function of the random vectors (Wi, . . . , W n ), of the random 
vector (TT51) and of their sum the following equality holds: 

0w+si(A)=0w(A)0s 1 (A). (34) 

We remember that the characteristic function of a random vector [X\, . . . ,X n ) whose joint distribution 
is a multivariate Gaussian 

^|exp(-i(Bx,x)) 

is (f>(\) — exp(— 1/2(B _1 A, A)) where B -1 is the covariance matrix of the vector (Xi, . . . , X n ). Thus the 

equality (jMf allows to determine the covariance matrix of the vector S 1 (cr) taking off the matrix A -1 
— i 

from 9~C G . By calculus it is easy to verify that 

Ok^ - A" 1 - x- (35) 

To complete the proof we have to show that the matrix x is positive definite. Consider the strictly 
convex function 

n , ,n \\ J ™ 

...,x n ) = In cosh ^a s Ji s x s + hi I ) = - ^ ona s Ji s xix s - G(x 1} ...,x n ) (36) 



2 

1=1 v s=l l.s=l 
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we can write IK^ = D a D Q JD a D a — where is the Hessian matrix of the function $. Since 
A = D a JD a we get "Kg = A — where 5€$ = D~ 1 JK$D a 1 . Multiplying x by the positive definite 
matrix we obtain 

Mc(A*)x = (A - ((a - 1 - A" 1 ) = MftOiJA" 1 . (37) 

Since all matrices involved in (J3TJ) are symmetric and JK$(/x)A~ 1 is positive define it follows that also x 
is positive definite. Hence the random vector (|18j). converges to a multivariate Gaussian which covariance 
matrix is X- 

2.2 Proof of Theorem H 

In order to prove theorem [2] we introduce the conditional joint distribution (see [9]) of a configuration 
<r conditioned on the event (mi(tr), . . . , m n (cr)) £ i?(A*, d) 

fW,J,h,d{o-} = Zjy (j h d ) ex P (277^ Js ' s ^ + ( h ' s )) 1 s(M,d)(^)^s(s) (38) 
where l_B(/x.d) is the indicator function of the ball B(fi,d), vg denotes the distribution of the random 
vector (S'i(cr), . . . , S n (cr)) on (M. N , YliLi p(°~i)) an d Zn(3, h, d) is the normalizing constant. 
To prove the theorem we need the following lemmas. 

Lemma 4. Let fi be a global minimum point of the function G given by Let k be the homogeneous 

type of (i. Define f — min{G(x)|x £ M™}. Then there exists a positive number 5^ such that for any 
t £ R™, any k £ N, any S £ (0, 5^] and any bounded continuous function cf> : R™ — > R 

lim e-<*> Nl/a V) ( f[ Nl ) ^e^ f e -^(*)+<t,Ni/»x )0(x)dx 
1=1 "> B{n,5) 

= 4>( fJ ,)J^eiq ? (-G^-^—) + (t,x))dx (39) 

See appendix for the proof. 

Lemma 5. Let P/v,j,h,d{f} be the joint distribution of cr — (a\, . . . , ct/v). Let Vj be the random vector 
Wi W n \ , {Si(a)-N im S n (<r) - N nf i n \ 



where W ~ iV(0, A -1 ) and A = D a JD a is a positive definite matrix. 
Then 

e -<t,N» /" ex p( -^(Jx,x> + (t,NTx> )Jjv(x,/i,d)(ix 

J exp f - y (Jx,x) + (t,N' 1 x)^I N (x,fj,,d)dx 
where (•) BG d denotes the expectation value with respect to the conditional distribution 138\) and 

L N (x,n,d)= / exp((Jax, s) + (h, s))dus(s). (42) 

7{AeB(M,d)} 
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See appendix for the proof. 

We give the proof of theorem [5] for /z = fx 1 . The other global minimum points are handled identically. 
Fix t e R" we choose the number 6 P > 0, p = 1, . . . , P according to lemma[4] For all 5 € (0, S p ] 



lim e-^ 1/2k ^ P Uf\ Nl ) 1/2kp e^ [ e -NG^ + (t,^^) dx 

= l eXP (" G ^(^) +(t ' X> )" X ^ 



(44) 



Define 

J^-(x,/i,d)= / exp((Jax, s) + (h, s))dv s (s). 

It is easy to verify that: 

I N (x,fi,d)+I c N (x,n,d) = cxp(7V$( x )). (45) 
where the function $ is given by (|36|) . For any p = 1, . . . , P, any < 9 < 5 and any < 5 < 6 P , define 

K N (t,6,^,e,k p ) = (f[Ni) 1/2kp exp(iV/ - {t.NVakp^P)) 



i=i 

x 



/ exp(-^(Jx ! x> + (t,N 1 / 2 ^x>)^(x,/xP,0) ( ix (46) 



Since K" = S^ 1 , 5i) US^/i 1 , 5i) and ^(x, /x 1 , cf) = e^*) - 7^(x, /Lt\ d) by lemma (with 7 = l/2fci) 
after have multiplied numerator and denominator by (IIILi Ni) 1 ^ 2kl e N ' $ we obtain 

/ (t,v V2fcl >\ = ^(Mi./j 1 ,/;!) - KN^Sujj^d, fci) + MivCMi,^ 1 ,^, fci) 
Gd L w (0 > <yi 1 M 1 > *i)-^w(0,*i > M 1 ,d,*i) + -WAr(0,<yi,M 1 > d > *i) 
where the random vector \i/2ki is defined by (I40p . 

MMi,/! 1 .**) = (n^) 1/2fel exp(iV/-(t,N 1 / 2 '= lM 1 )) 

x / exp (- — (Jx,x) + (t,N 1/2fcl x) + AT$(x)W 

™ x l/2fci 

IJJV,) exp(7V/-(t,N 1 ^ V i)) 
(=1 

/ exp(-7VG(x) + (t,N 1/2fcl x))dx (48) 



X 



and 



z=i 



x/ expf- — (Jx,x) + (t,N 1 / 2fel x))/ A r(x,/x 1 ,(i)dx. (49) 



To proceed we need the following: 
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Lemma 6. Let /J, 1 , . . . , fi p be minima points of the function G given by Let ki, . . . , k p be their 
homogeneous type. For any 9 > there exists e > such that for each p — 1, . . . , P 

K N {t,S p ,fj, p ,e,k p ) = 0(e~ Ne ) asN^oo (50) 

The proof of the lemma is in the appendix. By ([6]) with p — 1 and 6 — d there exists ei > such that 

K N (t,S 1 ,n\d,k 1 )=0(e- N ^). (51) 

Now we will prove that there exists also 62 > such that 

M N (t,5 1 , l i 1 ,d,k 1 )=0(e- Ne '). (52) 



Define the closet set 



Then 



V = W l -[jB{fiP,5 p ) (53) 



B c ( fJ , 1 ,S 1 )(lVu{jB( fJ ,P,5 p ). (54) 



D=2 



Thus we can write 

M N (t,5i,ij},d, fci) 

l/2fc 



(=1 

X 



/ \ 1/ jLK\ 

= (!]>) exp(AT/- (t.N 1 /** m X » 

!=1 

/ exp(-^(Jx,x) + (t,N 1 / 2fel x))//v(x, /it ,d)dx (55) 



Since d < <5 and H/x 1 -£t p || > £ for p = 2, . . . , P we have B(fi 1 ,Si) C B c (p p ,5-d) hence for each x e M" 
and p = 2, . . . , P 

I N (x, t i 1 1 d)<I c N (x lt iP 1 S~d). (56) 

Moreover by |45|) 

I N {x,n\d) <cxp(iV$(x)) (57) 
Using (f56|) in the integrals over B(/j, p , 8 P ), p = 2, . . . , P and ([57)1 in the integral over V we obtain 

/ " \l/2fei 

Mat^,^,/! 1 ^^!) = (n^j exptATf-^N^'V)) 



x 



/ exp(-7VG(x) + (t, N 1 / 2 ^))^ 
p 

ex P (<t, N 1 / 2 ^ - /x 1 ))) ^(t, /lx p , S-d,h) (58) 



p=2 
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Since the set V does not contain minima points, for any t € R ra there exists e > such that 

e Nf [ exp(-AG(x) + ( t , N x / 2fcl x))dx = 0(e~ Ne ) N -> oo. (59) 

Applying (j5"9")l to the first term of the right-hand side of and lemma [3] to each term of the sum, the 
result ((55J holds. 

By ([HI]) and ([52) we have 



> a "* L N (0,S U fj,\h) 



where in the last identity we use ((33]) . By the assumption on the random vector W the theorem [2J is 
proved. 

3 Examples 

We now analyze the case of two populations of the same cardinality. The Hamiltonian 

H N (m 1: m 2 ) = --^(juTnl + J 22 to 2 + 2J 12 m 1 m2 + 4/iimi + 4/i 2 to 2 ) (61) 

is a convex function of the magnetizations if the reduced interaction matrix J is positive definite, that is 
Jn > and Jn J22 — J12 > 0. A stationary point (/xi,/x 2 ) of the function G 



G{xi,x 2 ) = ^{j\\x\ + 2Ji 2 xxx 2 + J 22 z 2 ) - ~ln (cosh (-y-^i + ~^ x 2 + ^1)) 
- In (cosh (^Y x i + -y x 2 + ^2)) 



(62) 



2 

is a minimum point of homogeneous type k = 1 if : 



(JC G )iiGui,/z 2 ) > 
det:K G (/•*!, At2) > 



(63) 



where: 



and the determinant is: 



(JC G )u( M i,M2) - \(2Jn - J 2 n(l - Mi) - ^2(1 ~ mI)) (64) 



det( Ml) /xa) = ^ (4 - 2J U (1 - M ?) - 2J 22 (1 - fi) + det J(l - - ■ (65) 

For example if we consider the particular case in which the external field hi and /i 2 are equal to zero and 
the parameters Jn and J 22 are the same, the stationary point (0,0) verifies the conditions (f6"3"| if: 

< J n < 1 I K Jn < 2 

U ^ (66) 

- Jn < Ji 2 < Jn I Jn — 2 < J12 < 2 — Jn . 
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Thus for these choices of the parameters the random vector 

'S'i(er) S 2 (cr 



(67) 



weakly converges to a bivariate normal distribution. To have a minimum point of homogeneous type 
k > 1 the Hessian matrix of G computed in minimum point ((11,112) must be equal to the matrix with 
zero elements. This condition means: 

Ju > 2 
J22 > 2 



J12 



-- 

Ju-2 



(68) 



J11 

2 _ J22 - 2 



fj- 2 = 



J22 



Only if the third partial derivatives of G computed in ((11,(12) arc equal to zero the point can be a 
minimum point. This is verified if and only if Ju = J22 = 2. Hence (111,112) — (0,0). Computing the 
partial derivatives of fourth order we can assert that (0, 0) is a minimum point of homogeneous type 
k = 2. Thus the limiting distribution of the random vector 

'Si(<r) S 2 (cry 



N: 



3/4 



N. 



3/4 



(69) 



"1 J, 2 

is proportional to e - 1 / 24 ( x i+xl) _ 

As the parameter J12 = 0, we are also able to describe the limiting distribution of the random vector 
(Si(a), S2(cr)), properly normalized, beyond the homogeneity hypothesis on the minimum points of the 
function G. That is why in this case S\(<t) and 52(er) arc independent random variables and thus their 
joint distribution is the product of the marginal distributions. We observe that as Ju — 2 and hi = 0, 
for each value of the parameters J22 and h2, the determinant of the Hessian matrix of G computed in 
the minimum point (0,7x2) is equal to zero. The type of the minimum point (0,7x2) is not homogeneous 
whenever the second partial derivative of G with respect to X2 is different from zero. This condition is 
verifies as J 2 2 7^ 2 or h 2 7^ 0. For these choices of the parameters, the limiting distribution of the random 
vector 

1 Si(cr) S 2 ((t)\ 



JV?' 4 ' VN~2 

is given by the product of a Gaussian distribution with an exponential distribution of the forth order. 



4 Conclusions and Outlooks 

In this paper we have generalized to multi-species Curie- Weiss models the study of the normalized 
sums of spins and their limiting distributions. We worked under a condition of convexity of the reduced 
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interaction matrix which allows us to use the Ellis- Newman method. The theorems presented in this work 
obtain a complete classification of the distribution when the first non vanishing partial derivatives are 
all the same order (homogeneity hypothesis). The extension to non convex interactions or the complete 
classification of the limiting distribution beyond the homogeneity hypothesis will be subject of further 
investigation. 

A Proof of lemma [2] 

Considering the function 

^ n n n 

G(xi,. ..,x n ) = - ^2 a t a s Ji s xix s - ^ aiU ( a s J is x s + hi] (71) 

l,s = l 1 = 1 s=l 

where x = (xi, . . . , x n ) G R" and t = (ti, . . . ,t n ) G { — 1, +1}", the following inequality holds: 

G(x,t) > minG(x,t). (72) 

Since the reduced interaction matrix J is positive definite the function on the right-hand side of (1721) is a 
quadratic function with a positive definite Hessian matrix. Thus by inequality (|72p the global minimum 
points (at least one) of the function G belong to a same compact level set. Moreover these points, 
solutions of the mean field equations (|12j) . are isolated and the function G is analytic. Then the number 
of global minimum points is finite. 
By (J72]) we have: 

[ exp(-NG(x))dx< [ exp ( - iVminG(x,t)W (73) 

JR" JR" V t / 

Since the argument of the integral on the right hand side of inequality (|73|) is a Gaussian density function 
the statement (|24|) follows. 

B Proof of lemma [3] 

Given 8±, . . . ,6 n real 



„, Wx S^-Nm W n S n (<r)-N n m n 

\( Nl y/2-y + (iVl) i- 7 - i<---' (iN yi/ 2 - T + (JV n )i-7 " 



P {v/^i^i + Si(tr) G Ei, . . . , JKw n + S n {cr) G E n } 



where Ei = (—00, (A^;) 1 ~<6i + iV/mj]. The distribution of the random vector (y/Ni Wi,..., v^Vn W n ) is 

<l,,,A ^ expf-i(Ax,x>) (74) 



(2tt)™ J * \ 2 
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where it is easy to verify that A = 1/NJ. We claim that since the matrix A is positive definite also A 
has this property. The joint distribution of the random vector (Si(er), . . . , S n (cr)) is: 

' exp ( — (Js,s) + (h,s) )dv s (s) (75) 



Z N (3,h) r \2N 

where ^s(s) is the distribution of (Si(cr), . . . , S n (cr)) on (R N , YliLi p( a i))- The distribution of the random 
vector PB]) is given by the convolution of the distribution ([71|) with the distribution (|75|) . Thus: 



P{ y/Ph W x + Si(ff) G -£■!> • • • i W„ + S n (rr) G E n } 



1 / det A 



Z N (3,h) V(2tt) 



x //n ex P (^( ~ (J(w-s),(w - s)) + (Js,s) ) + (h,s) )dvs{s)dx 



Ei xl" 
i=i 



1 /detA s 
Z N (3,h) {j2n^ J 

exp ^- ^( Jw > w )^) J ex P ^( Jw ' s ) + (h,s)^dj/s(s)dw 



i=i 



where Ei denotes the Cartesian product of the sets Ei . 
i=i 

Since 



jf exp ^(Jw,s) + (h,s)Jdi/5(s) 

= n / N exp ( n Gi ( ^ + M zZ j ip w p) ) n 
= n n / exp i^f ^ + j; p w p ) Vf^) 



making the following change of variables 



wi - Nimi , 



and integrating over s, we obtain: 
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P{ t/WI W x + S x {a) G E x , . . . , W n + S n (a) e £„} 
ft W 7 



/=i 



/ det A 
Z N (J,h) \J2^ 



J\-J exp (- y £ + m,) + 

n n v 

A^; In ( cosh (hi^ct p Ji. p (^-^y +m p ^^jdx 1 ...dx T 



1=1 p=i 



1-7 



i=l 



/det A\ 2 



Zjv(J,h) V(2tt)", 

x / . . . I exp (- NG(^j + mi,. . . +m n ^dx x . . .dx n . (76) 



— oo 



Taking X — >■ oo, . . . , 9 n — > oo the (|75|) gives an equation for Zjv(J, h) which when substituted back yields 
the result (j2"T)l . The integral in the last expression is finite by (|24l) . 

C Proof of lemma [4] 

To easy the notation we define 7 = l/2fc. Making the change of variable 

xi = m + j4y 1 = 1, ...,n (77) 

the left-hand side of (1391) becomes 



^i^ll</^ + ^ ^ { NB ^ + ^''^ + (t ' X> )" X (78) 
where B is the function defined in (l29l) . By the conditions expressed in (|30|) and dominate convergence 
theorem the limit ([75)1 is equal to 



t,u) )du (79) 



Since G^. ( * 2k ) > for every x different from zero, the integral in (JTOJ) is finite. This completes the 
proof of the lemma. 
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D Proof of lemma [5] 



By following the same proof of lemma |3] we have that the distribution of the random vector V~ ( 
given by 

exp (- § (J(> + #f)) / w(M+ m,t*,d)dx 



N x 



N 7^"~ N 



Sf 7 V 



where 



Thus 



= / exp((Ja(/x + -^-),s) + (h,s))^ s (s). 

V JN7 7 '{#eB(M)} JN7 



exp ( - ir(J(M + ^r), ^ + =?rW 



<t,x> , 



Making the change of variable 



the statement (|4T1) holds. 



ui = m + j^y l = l,...,n 



E Proof of lemma [6] 



We prove the theorem for p = 1. The proofs for other p are similar. We observe that 



iVi 



;=i 



u 

i=l 



Sl(<r) „i 



AT; 



> 



AT, 



>n} + e 



where 9 = y/nQ. Thus 



I c N (x,fi\9) < V ( ! exp((Jc*x,s) + (h,B))du s (fi) 

i =1 \->{-%>-ri+S} 



exp((Jax, s) + (h, s))<fz/g(s) 



Consider one of the integrals of expression (|85|) 

exp((Jo;x, s) + (h, s))dfs(s) 



L 



;>»]+(>} 



Y\ COSh y ^ UqJlqXq + h 
1=2 q=l 



Ni 



/ exp(5i(<r)( ^a q J lq x q + hi\ ) TT dp(oi). 

'{Si( -)>Jv 1 ( M ;+e)} v //xx 



9=1 
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By Chebishev's inequality for any r > 



„ n 

/ exp[S 1 (cr)( ^S^aqJiqXq + hi)) T\ dp(<Ji) 

J {s 1 ( ( t)>n 1 (^ 1 +S)} v v ^ ieA 



< exp(-aiJiir7Vi(/z} + 0)) 

x / exp(ai Jut ^ cri) exp(^ cr^ai Jn^i + a 2 Ji2^2 + 

= exp(7Vi(-ai Jnr(^J + 0) + lncosh(ai Jn(xi + r) + a 2 Ji2X2 + hi))). (86) 



By the mean field equations (|T2|) we have: 

- — (lncosh(ai Jnxi + a 2 Ji2^2 + /ii))^ 1 ) = ceiJufa. (87) 

Thus we can choose 8 > and r > sufficiently small such that 6 < Si and 

lncosh(o!i Jn(%i + t) + 02^12^2 + hi) < lncosh(ai JuXi + c*2 J12X2 + hi) + cti Juh\t + —oliJutO (88) 

for each x 6 B{fj},8). The other integrals in (|85p are handled in a similar way. At last applying the 
bounds (ESI) and flSHJl to ([55|). for all x e B(/x\<5) we obtain: 

I c N (x,fj,\9) < 2nexp (jv( - ^a 2 J U T0 + $(x))) (89) 

where $ is given by (|56"]) . Hence 

^(t^juSMi) < 2n(7V ; ) 1/2fel exp (iV/ - (t, N 1 / 2 ^ 1 ) - y«?Jiirff) 

x / exp(-iVG(x) + (t, N x / 2,!l x))dx. (90) 

Applying the lemma H] to the expression on the right-hand side of (|90| we obtain: 

K N (t,S,n\9,ki) = 0(e-% a i JllTS ) asiV^oo. (91) 

We now bound 

K N (t, Sun 1 , 9, hi) - K N (t, S, /x\ 9, k x ) 

= (iV ; ) 1/2fcl exp(7V/ - (t,NV^ M i)) 

x / exp(-^(Jx ! x) + (t,N 1 / 2fel x))^(x,/x 1 ,^x. (92) 

•/B(/*Vi)--B(mM) V 2 / 

By (j45|) we have: 

I C N (*,»\0) < exp(7V$(x)). (93) 

By definition of the function $ we get 

Jr*(t, Sijfi 1 ^, h) - K N (t, S, n\9, hi) < (iVi )i/2 ^/^t.NV-V) 

x f exp(-iVG(x) + (t,N 1 / 2fcl x))dx. (94) 
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Making the change of variable 



we obtain 



Xl 



N, 



where 



Observing that 



E = 



{ 



E C 



<*}\{ 



u 



Nl/2fci 



N!/2fel 

< 4 }' 



< 5 



it follows that as N —> oo, for some eo > 

K N (t, 6 uf i\6, ki) - K N (t, 6, fi\e,h) = 0{e m ") 
The statement ^SDJ) follows by (|9"Tj) and 



(95) 



JGv(t, gi, /i 1 , 0, fci) - K N (t, S, n\ 9, h) < jf exp f - i + ^i/ 2fcl ? ^)) + u )) du ( 96 ) 



(97) 
(98) 

(99) 
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